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Abstract. 

We present a method for computing the partition function of a caloron ensemble taking into account the interaction of 
calorons. We focus on caloron-Dirac string interaction and show that the metric that Diakonov and Petrov offered works well 
in the limit where this interaction occurs. We suggest computing the correlation function of two polyakov loops hy applying 
Ewald’s method. 

Keywords: QCD, calorons, interacting ensemble, Ewald’s method, confinement 
PACS: ll.15.Ha, 12.38.Aw, 12.38.Lg, 12.39.Pn 


INTRODUCTION 


KvBLL instantons or calorons with non-trivial holonomy can be taken as the degrees of freedom in Yang-Mills vacuum 
in any temperature. The non-interacting ensemble of calorons has been studied by Diakonov and Petrov lH. They 
obtained the heavy quark potential by calculating the correlation function of two Polyakov loops: 

= (TrL{r)TrL^{r)). (1) 

To compute the correlation function, the partition function of caloron ensemble has to be calculated. Therefore, one 
has to hnd the contribution of a caloron in the partition function. Diakonov and Gromov iQl computed this contribution 
using the metric of the caloron moduli space in terms of N different SU(N) BPS monopoles 

Z = e-^^' j Da^DxDx^xp j d^a^W^^v- j • 


Where Yp represents the caloron collective coordinates and the Jacobian J is the determinant of the moduli space metric 
tensor. Diakonov and et a/.Q computed the full partition function for KvBLL instanton in SU(2) gauge theory: 


Z = f d^zid^Z 2 dffT^c(^^y exp [-VP (v) - iTZrnP'' (v)] 
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where C is a constant and 
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v^v^P"(v) = - V][v- 7r7'(l - —)]. 






Using the partition function Z, for computing the correlation function of Polyakov loops in eq.[Tl one can obtain the 
potential between a quark-anti quark pair. 

In this research we try to add the contribution of the interacting calorons to the calculations. As a result of this change, 
the metric of the interacting calorons moduli space should be modihed. However, one may keep the same moduli 
space with the same metric as the non interacting calorons, but introducing a potential on the moduli space which 
corresponds to the interaction between calorons. The main task of this research is to give some ideas about computing 
the partition function of the interacting calorons and hnally calculating the potential between static quarks. However, 
one can compute correlation function of the Polyakov loops without using metric by applying Ewald’s method. 

In the next two sections we use the sum ansatz idea and the caloron-Dirac string interaction to suggest how one may 
calculate the partition function. In section 4, we discuss about the caloron-anticaloron ansatz and in section 5, we argue 
about the interval where the metric of Diakonov and Petrov is positive dehnite and we explain that our calculations in 
section 3 is done in this interval. In section 6, we apply Ewald’s method to caloron ensemble for large caloron radius 
limit. 




SUM ANSATZ 


Even though, two calorons do not interact at large distances in algebraic gauge but they may interact when they get 
close to each other. 

Consider two calorons with non-trivial holonomies. Because A 4 ’s of these two calorons do not vanish even at spatial 
infinity, calorons correlate with each other at infinite separation. Therefore one should not simply add two single 
caloron gauge helds. It is possible to choose algebraic gauge, in which all components of potential of a caloron vanish 
at spatial inhnity, and then add the potentials of calorons: 

(3) 

We recall that to compute the partition function of sum ansatz, A^’s in eq. [3 should be transformed back to the 
periodic gauge. This gauge transformation is dependent to the holonomy, therefore two calorons should have identical 
holonomies. The held strength of this ansatz can be written as 

(a(') +a(2)) + V (aP)) ^ 

where the additional term (the third one on the right) is due to the non-linear nature of the theory. Calculating the 
action, we have an extra term in addition to the action of individual calorons that can be interpreted as the effective 
classical potential corresponding to the interaction between two calorons 


exp(-5;„,) = exp 
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Now we should write the partition function of two interacting calorons. We should note that in computing the partition 
function we change the dominant path to combination of two calorons in periodic gauge, so we must consider the 
quantum fluctuation about this new path and compute the metric and quantum potential. But as the hrst order of 
approximation we can use the metric of one individual caloron and write the partition function as: 
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where Z, has been derived in eq.|2]. 

Now we study the ensemble of K calorons. To simplify the calculations we restrict ourselves to a system which consists 
of many pairs of calorons. The calorons in each pair interact with each other but each pair does not interact with the 
others. The grand canonical partition function of these K pairs can be written as: 

K K 

Where f is the fugacity and should be calculated, and Z 2 c can be obtained from eq.|5] With this partition function, one 
can calculate the correlation function of Polyakov loops and then compute the heavy quarks potential. 

In the next section, we argue about another alternative to include the interaction between calorons in the partition 
function. 


CALORON-DIRAC STRING INTERACTION 

The next idea for computing the partition function of interacting calorons is using the caloron-Dirac string interaction 
which has been studied by Gerhold et al. |3l. The interaction between the caloron and the Dirac string had not been 
taken into account in the non interacting caloron partition function and therefore this is a new effect in computing 
the partition function of interacting calorons where two calorons are located on top of each other. It means that two 
calorons interact with each other by their Dirac strings. 

In the large caloron radius limit, p ^ j3 , where p is the radius of the caloron, the vector potential of the caloron 



in algebraic gauge is dominantly Abelian and governed by the third component in color space. For convenience, the 
monopole axis is taken along the third axis. 



(x) 


<Ag(x), 





This dominant component is called a Dirac string. 

Since the interaction due to the Dirac string is not negligible, when an object (monopole or caloron) is placed on the 
top of the Dirac strin^f a caloron, the ordinary sum ansatz can not be used for the superposition of caloron and this 
object. Gerhold et a/.|3] derived the superposition by removing the Dirac string in the algebraic gauge doing a proper 
gauge transformation (G), then they added the gauge fields of this caloronand that object and finally applied the inverse 
gauge transformation to have the original Dirac string back. The result is; 

(x) = (x) + (x) ^ (6) 

Therefore the effect of the Dirac string is included by adding an unchanged caloron to a gauge rotated object. 

Now consider two interwined calorons such that one of the monopoles of a caloron is located on the top of the Dirac 
string of the other caloron and vice versa. In this case one should rotate only the monopole of each caloron which 
touches the Dirac string of the other caloron as explained in eq. |6] These rotations do not affect the Dirac strings, 
because monopoles rotate around the monopole axis which is along the third axis. We should recall that the Dirac 
strings are also defined along the third axis. Applying these rotations, one can repeat the approach of the previous 
section and obtain the field strength of and use eq.|4]to compute the interaction term. Again, we restrict ourselves 
to a system which consists of many pairs of calorons. The calorons in each pair interact with each other but each pair 
does not interact with the others. We can compute the partition function of K non-interacting pairs of calorons as in 
previous section. 


CALORON-ANTICALORON ANSATZ 

The sum ansatz idea can be applied to a caloron and anticaloron system. Although the combination of solutions of 
caloron and anticaloron is not the solution of equation of motion but both of these objects minimize the action and 
may have identical integral measure and partition functions. They may be saddle points of the partition function. In 
1984, Diakonov and Petrov lUt] argued about the instanton-anti instanton solutions and wrote:"Thus, from physical 
considerations it is clear that the main contribution to the partition function should be given not by exact classical 
solutions, but by approximate solutions of the instanton-anti instanton type. This should not surprise us since the 
approximate solutions may have a larger entropy (= statistical weight) than the exact solutions which correspond to the 
local minima of the action. Anyhow, one is interested in maximizing the partition function of a theory, and instanton- 
anti instanton configurations may well give a larger contribution to it than those of multi-instantons." Therefore, there 
is a possibility of constructing the partition function from caloron anticaloron pairs instead of calorons or anticalorons 
only. However, we recall that the caloron and anticaloron should have the same holonomy if one wants to do an 
appropriate gauge transformation, from algebraic gauge to the periodic gauge, on the combination.. 

Since the function G which is applied to remove the Dirac string of anticaloron is the same as the caloron’s |0], we 
can apply the caloron-Dirac string interaction to the caloron-anticaloron ansatz. Thus, we can consider an ensemble 
with equal numbers of calorons and anticalorons in which one caloron interacts with only one anticaloron (either 
interaction due to the ordinary sum ansatz or caloron-Dirac string interaction). Then, we again assume that these pairs 
do not interact and we use the same procedure as the previous section to calculate the partition function. 


THE METRIC OF MODULI SPACE 


In 2009, Bruckmann et al. |@] showed that the metric of moduli space of dyons which Diakonov and Petrov |01 offered, 
does not have the essential requirement of positive definiteness, throughout the whole configuration space. They argue 
that this metric is positive definite only for the low dyon density. 

We recall the approximate metric of the same kind dyons in SU(2) is 10]: 


G = 








The eigenvalues of this metric are: 
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where d is the dyon separation. Det (G), which represents the weight factor, is positive for d > ^ , where T shows 

the temperature. Therefore for the dyon separation smaller than this critical distance, where two dyons overlap, this 
approximate metric does not work and one should refine it. (The metric of different kind dyons is positive definite for 
all dyon separations, because the minus sign in the above formula changes to plus sign and therefore the eigenvalues 
of the metric are positive for all d). But we recall that for the caloron-Dirac string interaction, we take the caloron in 
the large raduis limit, p ^ j3 , where the vector potential of the caloron is dominantly abelian. Rewriting the limit of 
positive definiteness. 



yielding ^ ^ ~ 0-45 . Therefore in the limit where we need to consider the Dirac string, the metric which Diakonov 

and Petrov offered is positive definite and works well. We shoule notice that for calculating the partition function from 
sum ansatz, one should either consider the calorons in the large raduis limit or modify the metric. 


EWALD’S METHOD 


To deal with the problem of positive definiteness, Bruckmann et al. isj] used the Ewald’s method to compute the 
correlation function of two polyakov loops in eq[T] for non-interacting ensemble. The vector potential for 2K dyons 
ensemble in the abelian limit is 
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where = —qri^ydvln{r — z). 

Then, one rewrites the polyakov loop 
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where <I>(r)= ^ ^ is a long-rang function. For maximal non-trivialholonomy, both dyons of the SU(2) gauge 

i=lm=l ' ' ' 

group have the same action. Therefore 
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To use Ewald’s method, one should mimic the infinite space by sampling the system to the so-called "super cell" of 
volume 1 } with finite dyon density. Each super cell is represented by n S , therefore <t>(r) can be written as 




neZ" ^ 


— Tj — nL\ ’ 


where j = (i, m). Then, one should split the above term to a "long-range part" and a "short-range part". The short-range 
part vanishes and the Fourier transform of the long-range part converges at spatial infinity. Back to our work, we have 
calorons instead of dyons in the limit of p ^ j3 . The vector potential of the caloron in this limit is 


A^(r) = --(J3fj3^(9vln‘P(r), 






where 'l'(r) = and zi and Z 2 represent the positions of two constituent monopoles. If the observing 
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point is far from the caloron, then r ^ Zi,Z 2 - In the first approximation, where ^ and ^ {i = 1,2) are negligible, we 
can write the vector potential as 
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Since we have a long-range vector potential, we can use the Ewald’s method. 

Now we can consider the two interacting calorons as one object and write its vector potential as in eql6]. This vector 
potential is also long-range and we can use Ewald’s method for an ensemble of K non-interacting two-caloron objects 
with appropriate separation and repeat the procedure of this section two compute the correlation function of Polyakov 
loops. 


CONCLUSION 

In this research we try to suggest computing the potential between interacting calorons with non-trivial holonomies.The 
caloron-Dirac string interaction may be used to obtain the interacting potential between two calorons. Having this 
potential in hand, we have two ways to compute the correlation function and then the heavy quark potential. One way 
is calculating the partition function by computing the metric of the moduli space of the caloron ensemble and the other 
one is using the Ewald’s method to include the effect of the long-range vector potential of calorons in the limit where 

p>j3 . 
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